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Exact solutions are found to the Dirac equation for a combination of 
Lorentz scalar and vector Coulombic potentials with additional non- 
Coulombic parts. An appropriate linear combination of Lorentz scalar 
and vector non-Coulombic potentials, with the scalar part dominat- 
ing, can be chosen to give exact analytic Dirac wave functions. The 
method works for the ground state or for the lowest orbital state with 
I = j - |, for any j. 

In a previous letter |l|], simple exact solutions were found for the Dirac 
equation for the combination of a Lorentz vector Coulomb potential with 
a linear confining potential that was a particular combination of Lorentz 
scalar and vector parts. In this letter, we extend the method of Ref.[|lJ to 
the more general case of an arbitrary combination of Lorentz scalar and 
vector Coulombic potentials with a particular combination of Lorentz scalar 
and vector non-Coulombic potentials, S(r) and V(r). The non-Coulombic 
potentials can have arbitrary radial dependence, but must be related by 



where E is the bound state energy for a particle of mass m. This is the 
same relation between the vector and scalar non-Coulombic potentials as 
was required in Ref.JlJ 




(1) 



m 
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The Dirac equation we solve is 



a ■ p + (3m - iA±M + V (r) + /3S(r) 



ip = Ei/j, 



(2) 



where cx and (3 are the usual Dirac matrices. The four component wave 
function ip can be written in terms of two component spinors u and v as 



ip = N 



a 



(3) 



with N an appropriate normalization constant. The two component spinors 
u and v satisfy the equations 



(a ■ p)v 



[a ■ p)u 



E-m + (A + 7/) - V{r) - S{r) 



E + m + 



r 

(A -77) 



V(r) + S{r) 



ii 



v. 



(4) 
(5) 



The key step in generating relatively simple exact solutions of the Dirac 
equation is to choose a particularly simple form for the function v(r) 



v{r) = 27(17 • r)u(r), 



(6) 



where 7 is a constant factor, to be determined by the solution to the Dirac 
equation. This is the form of v(r) that was found in Ref.|l[ for a Coulomb 
plus linear confining potential. This ansatz for v(r) has also been used as 
the basis for generating approximate saddle point solutions for the Dirac [0 
and Breit|| equations. 

For now, we limit our discussion to spherically symmetric states u(r) that 
depend only on the radial coordinate. We will discuss orbitally excited states 
later in this paper. Using the form of v(r) given by Eq. (|B|), equations (4) 
and (5) reduce to two first order ordinary diffential equations for u(r) 



du 1 

dr 7 
du 
dr 



E-m- V{r) - S{r) + 



(A + 77 — 2 7 ) 
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E + m - V(r) + S(r) + 



r 

(A - V ) 



u(r) 



u(r). 



(7) 
(8) 
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Equations (7) and (8) are two independent equations for the same quantity, 
so that each term in one equation can be equated with the corresponding 
term in the other equation having the same radial dependence. This leads to 

2 _ m - E _ S(r) + V{r) _ X + r] - 27 
7 ~ m + E ~ S{r) - V{r) ~ 77 — A ' ' ' 

The relations in Eq. (9) can be rearranged, after some algebra, to give 

7 = ££ (10) 

with 



b = ±^l-\ 2 + r] 2 . (11) 

The constant b can have either sign. Although h must be positive in the 
pure Coulombic case, we will see that a negative b is possible if the Lorentz 
scalar potential S(r) is more singular at the origin than 1/r. The bound 
state energy can be written as 

or as 

The wave function u(r) can be found by solving differential equation (8) 
to give 

(14) 



u(r) = r b 1 exp 



-a H J S{r)dr 



where the constant a is given by 

a = 1 {m + E) = ±Vm 2 - E 2 . (15) 

The constants 7 and a could be negative if S(r) approaches a large enough 
negative constant or diverges negatively as r becomes infinite. The integral 
in equation (|l4l) can diverge at the origin or as r — > 00, or for any finite r, 
as long as the quantity in square brackets in Eq. ([14]) remains negative. 

The constraint equation (|i~3|) shows that in order for this class of exact 
solutions to apply, the Lorentz vector and scalar non-Coulombic potentials 
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must have the same radial dependence. As long as this constraint is satisfied, 
the results in Eqs. (10)-(15) represent a complete exact solution for the wave 
function and bound state energy of the Dirac Hamiltonian given in equation 
(2). 

We note from Eq. (|T2|) that the energy seems to depend only on the 
Coulombic coupling constants A and rj, and does not seem to depend on the 
non-Coulombic potentials. However, we could alternatively say from Eq. (|13|) 
that the energy depends only on the ratio of the non-Coulombic potentials 
V(r) and S(r), and does not seem to depend on the Coulombic coupling 
constants. Then, Eq. ( |T2| ) could be considered a constraint equation on the 
Coulombic coupling constants. The actual situation is that the energy does 
depend on all the potentials but, because of the severe constraints imposed 
by equations ( |T2D and (|13|) taken together, the energy can be written in terms 
of one set of potentials or the other. Although the possibility of this class 
of exact solutions is limited by the constraints on the potentials, this still 
permits a wide range of non-Coulombic potentials. 

We now consider conditions imposed on the potentials and the wave func- 
tion parameters by the physical requirements that the potentials be real and 
the wave function normalizable. We see from Eq. (|9]) that 7 must be real, 
and then from Eq. (|10"D that b must be real. This requires the Coulombic 
potentials to satisfy the condition 

l-A 2 + r? 2 >0. (16) 

The reality of 7 restricts possible bound state energies to the range 

— m < E < m. (17) 

Note that negative energies can occur, but E+m cannot be negative. Also, E 
cannot equal ±m, because this would lead to an unormalizable wave function. 
This condition on E, along with Eq. flT3|), means that V(r) must always be 
less in magnitude than S(r). 

We discuss the remaining conditions on the parameters in terms of three 
sub-classes of solution: 

1. The "normal" class of solutions has b, 7, and a all positive. In this 



case, we see from Eq. (10) that the Coulombic potentials must satisfy 
the further condition 

A + 77 > 0. (18) 
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2. Sub-case 2 has b negative, with 7 and a still positive. The constant b 
can be negative if the product aS(r) is positively divergent at the origin 
faster than 1/r. Then each of Eqs. (10)-(15) holds just as for positive 
b, and the wave function is still normalizable. Sub-case 1 with positive 
b transforms smoothly into the pure Coulombic solution as the non- 
Coulombic potential tends to zero everywhere. But this is not true for 
sub-case 2 with b negative. This sub-case requires the non-Coulombic 
potential to be dominant at the origin, and so has no corresponding 
pure Coulombic limit. 

3. Sub-case 3 has a negative 7 and a negative a, while b can have either 
sign, as discussed in sub-cases 1 and 2 above. A negative a is possible 
if the non-Coulombic potential diverges or approaches a constant as 
r — > 00, so that the integral in Eq. (|I5|) diverges faster than r at large 
r. Since a is negative, the potential S(r) must be negative at large 
r. Then all of Eqs. (10)-(15) hold as for positive a, and the wave 
function is still normalizable. This case is highly unusual, because it 
allows the possibility of a potential that is negative everywhere and 
diverges negatively at both the origin and infinite r. We know of no 
other example in quantum mechanics where a potential that diverges 
negatively at infinity can lead to a normalizable ground state. The 
reason this is possible here can be seen from Eqs. (7) and (8). There 
it is seen that S(r) enters the differential equations for u(r) only in 
the combinations 76" (r) or S(r)/j. Since these effective potentials are 
positive, the resulting wave function is normalizable. As with sub-case 
2, the case with 7 and a negative cannot not approach a pure Coulombic 
case because the non-Coulombic potential must be dominant at large 
r. 

We now look at some special cases. If the non-Coulombic potentials are 
absent, then the solutions are for a general linear combination of Lorentz 
vector and Lorentz scalar Coulombic potentials. If either constant, A or 77, is 
zero, we recover the usual solutions of the Dirac equation for a pure scalar 
or vector Coulombic potential. The Coulombic potentials cannot both be 
absent (while keeping a non-Coulombic part) because then 7 would be zero 
and b one, leading to a constant, unnormalizable wave fucntion. 
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For a power law non-Coulombic potential of the form 
S(r) = fi(s + l)r s , 
the wave function will be given by 



u[r\ 



r b l exp 



- a I r — —r s+1 
m 



(19) 



(20) 



All three sub-classes are possible for this wave function, depending on the 
ranges of the parameters b, a, /i, and s. 

The method does not work for radially excited states, because the simple 
ansatz of Eq. (Q) for v(r) does not lead to consistent equations for du/dr in 
that case. But, the method does work for the lowest orbitally excited state 
for which I — j — |. This can be seen by writing the Dirac wave function in 
terms of radial functions f(r) and g(r) 



ij = N 



f(rW i 

J .2 

-ig(r)y J i 



(21) 



where y\ is a two component angular spinor function corresponding to total 
angular momentum j and orbital momentum I. The radial functions satisfy 
the equations 



d_ 

dr r 

' d_ 1+- 

, dr r 



K 



f 



E + m + ^ — ^ + 3{r) - V(r) 



K 



E — m + 



r 

(A + 
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S{r) - V(r) 



(22) 
(23) 



where k = j + \ is the principal quantum number of the state. 
For the radial function g(r), we make the ansatz 

g(r) = 7/(r). 



(24) 



Then, equations (22) and (23) reduce to two first order differential equations 
for f(r) 



d£ 
dr 
d£ 
dr 



1 

? 



-7 E + m- V(r) + S{r) + 
E-m- V(r) - S(r) + 



(X-r,) (k-1) 
(A + 77) 7(« + l) 



f(r) (25) 
f(r). (26) 
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Equating the corresponding terms having the same radial dependence in 
these two equations results in 



2 m — E S(r) + V(r) A + rj — 2'jk 
m + E S(r) — V(r) rj — A 

Equation (|2"7|) is the same as Eq. (|S|) with the replacements 

A — > X/k, rj — ► 77/ 7c. 



(27) 



(28) 



Thus, equations (10)-(13), and (15) hold for the orbitally excited case, with 
the replacements A — > A/k and rj ^ rj/n. The radial wave function is /(r) is 
given by 



„k6-1 



—a (r H / S(r)dr 

\ m J 



(29) 



Note that the constant 6 used in this paper differs from the b defined by Eq. 
(22) of Ref. 1. The b in this paper is the b in Ref. 1 divided by k. 

In summary, we have presented a class of potentials for which the Dirac 
equation has relatively simple exact analytical solutions for the ground state 
of each angular momentum state with I — j — 5. 
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